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The computational complexity of Dezert-Smarandache Theory (DSmT) increases exponentially with the 
linear increment of element number in the discernment frame, and it limits the wide applications and 
development of DSmT. In order to efficiently reduce the computational complexity and remain high 
accuracy, a new Evidence Clustering DSmT Approximate Reasoning Method for two sources of information 
is proposed based on convex function analysis. This new method consists of three steps. First, the belief 
masses of focal elements in each evidence are clustered by the Evidence Clustering method. Second, the 
un-normalized approximate fusion results are obtained using the DSmT approximate convex function 
formula, which is acquired based on the mathematical analysis of Proportional Conflict Redistribution 
5 (PCRS) rule in DSmT. Finally, the normalization step is applied. The computational complexity of this 
new method increases linearly rather than exponentially with the linear growth of the elements. The 
simulations show that the approximate fusion results of the new method have higher Euclidean similarity 
to the exact fusion results of PCR5 based information fusion rule in DSmT framework (DSmT + PCR5), 
and it requires lower computational complexity as well than the existing approximate methods, especially 
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1. Introduction 


As a novel key technology with vigorous development, informa- 
tion fusion can integrate multiple-source incomplete information 
and reduce uncertainty of information which always has the con- 
tradiction and redundancy. Information fusion can improve rapid 
correct decision capacity of intelligent systems and has been suc- 
cessfully used in the military and economy fields, thus great atten- 
tion has been paid to its development and application by scholars 
in recent years [1-9]. As information environment becomes more 
and more complex, greater demands for efficient fusion of highly 
conflicting and uncertain evidence are being placed on information 
fusion. Belief function theory (also called evidence theory) referred 
by Dezert-Smarandache Theory (DSmT) [9] and Dempster-Shafer 
Theory (DST) [10,11] can well deal with the uncertain and conflict 
information. DSmT, jointly proposed by Dezert and Smarandache, 
is considered as the general extension of DST, since it beyonds 
the exclusiveness limitation of elements in DST. DSmT can ob- 
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tain precise results for dealing with complex fusion problems in 
which the conflict is high and the refinement of the frame is not 
accessible [9]. Recently, DSmT (belief function theory) has been 
successfully applied in many areas, such as, Map Reconstruction of 
Robot [12,13], Decision Making Support [14], Target Type Tracking 
[15,16], Image Processing [17], Sonar Imagery [18], Data Classifi- 
cation [19-21], Clustering [22,23], and so on. Particularly, the very 
recent credal classification methods [12,15,20] working with belief 
functions have been introduced by Liu, Dezert, et al. for dealing 
with uncertain data, and the object is allowed to belong to any 
singleton class and set of classes (called meta-class) with differ- 
ent belief masses. By doing this, the credal classifiers are able to 
well capture the uncertainty of classification and also efficiently 
reduce the errors. However, the main problem of the application 
(e.g. classification task) of DSmT is that when the focal elements’ 
number increases linearly, computational complexity increases ex- 
ponentially. 

Many approximate reasoning methods of evidence combination 
in DST framework were presented in [24-26]. But these methods 
cannot satisfy the small amount of computational complexity and 
less loss of information requirements at the same time. In recent 
years, there are some important articles [27-34] dealing with the 
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computational complexity of the combination algorithms formu- 
lated in DSmT framework in different ways. Djiknavorian [27] has 
proposed a novel method and a Matlab program to reduce the 
DSmT hybrid rule complexity. Smarandache has proposed Smaran- 
dache’s codification [35] to describe the parts of a Venn diagram in 
DSmT framework, which is easy to read and apply in unions and 
intersections of sets [28]. For manipulating the focal elements eas- 
ily, Martin [28] has proposed a Venn diagram codification, which 
is more complex for the readers than Smarandache’s codification. 
However, the DSmT (belief function theory) complexity can be re- 
duced by using this practical Martin’s codification and only consid- 
ering the reduced hyper-power set D° after integrating the con- 
straints in the codification at the beginning of the general belief 
function framework proposed in [28]. Abbas [29,30] has proposed 
a DSmT based combination scheme for multi-class classification 
which also reduces the number of focal elements. Li [31] has pro- 
posed a method for reducing the information fusion complexity, 
which is different from the above methods by reducing the com- 
bined sources numbers instead of reducing the number of focal 
elements. Li and other scholars [32-34] also proposed an approx- 
imate reasoning method for reducing the complexity of the Pro- 
portional Conflict Redistribution 5 (PCR5) based information fusion 
rule within DSmT framework. However, when processing highly 
conflict evidences by the method in [32], the belief assignments 
of correct main focal elements transfer to the other focal elements, 
which leads to low Euclidean similarity of the results in this case. 

Aiming at reducing the computational complexity of PCR5 
based information fusion rule within DSmT framework (DSmT + 
PCR5) and obtaining accurate results in any case, a new Evidence 
Clustering DSmT Approximate Reasoning Method for two sources 
of information is proposed in this paper. In Section 2, the basics 
knowledge on DST, DSmT and the dissimilarity measure method 
of multi evidences are introduced briefly. In Section 3, mathemat- 
ical analysis of PCR5 formula is conducted, which discovers every 
conflict mass product satisfies the properties of convex function. 
A new DSmT approximate convex function formula is proposed 
and error analysis of the proposed formula is also presented. Based 
on the error analysis, an Evidence Clustering method is proposed 
as the preprocessing step and the normalization method is applied 
as the final step of the proposed method for reducing the approx- 
imate error. The process of the proposed method is given, then 
analysis of computation complexity of DSmT + PCR5 and the pro- 
posed method are presented. In Section 4, the results of simulation 
show that the approximate fusion results of the method proposed 
in this paper have higher Euclidean similarity with the exact fu- 
sion results of DSmT + PCR5, and lower computational complexity 
than existing DSmT approximate reasoning method in [32]. The 
conclusions are given in Section 5. 


2. Basic knowledge 


In this section, we will give an overview of the basics knowl- 
edge on DST and DSmT, which are closely related to our work in 
this paper. 


2.1. Dempster-Shafer Theory (DST) 


Let us consider a discernment frame © = {6}, 62,---,@n} con- 
taining n elements 6;,62,---,), which is the refinement of the 
discernment based on Shafer’s model. The basic belief assignment 
(bba) is defined over the power-set 2© which consists of all sub- 
sets of ©. For example, if one has © = {6}, 62, 03}, the power set is 
given by 2° = {O, 01, 62, 63, 01 UO2, 61 U63, 62 U63, 6; U62 U63}, and 
the bba m(.) : 2° — [0, 1] on the power set is defined by [10,11] 


m(X;)=0, xXj=9 (1) 


> m(Xj) =1 (2) 


Xj€2©,1<i<n 


The element Xj; is called focal elements if it holds m(X;) > 0. 
Dempster’s rule is often used for the combination of multiple 
sources of evidence represented by bba’s in Shafer’s model, and 
it requires that the bba’s must be independent. The bba of the ith 
source of evidence is denoted m;. The Dempster combination rule 
is defined by Equation (3) and the conflict in Dempster combina- 
tion rule, denoted by C, is defined by Equation (4) [10,11] 


1 
mps(Z)=7——_ Dm (Xi)-ma(Xj) VCO (3) 
XjNX j=Z,iAj 
C= SY) mj (Xi) -m(X)) (4) 
Xi,XjCO,iAj 
XjNXj=0 


One can see that all the conflicting beliefs C has been redis- 
tributed to other focal elements. Dempster’s rule usually produces 
very unreasonable results in the fusion of high conflicting infor- 
mation due to the redistribution of conflicting beliefs. In order to 
solve this problem, many alternative combination rules like Pro- 
portional Conflict Redistribution 1-6 (PCR1-6) rules [36,38,39] have 
been developed. 


2.2, Dezert-Smarandache Theory (DSmT) 


DSmT [29] overcomes the exclusiveness limitation in Shafer’s 
model. In many fusion problems, the hypotheses can be vague in 
reality and the elements are not precisely separated which don’t 
satisfy Shafer’s model. The hyper-power set denoted by D® is built 
by applying operator M and U to the elements in © [36,37]. Let 
us consider a simple frame of discernment © = {6;, 62}, then one 
gets D© = {D, 0, 02, 63, 0 UO2, 016}. The bba in DSmT is defined 
over the hyper-power set as m(.) : D® — [0, 1]. 

In the combination of multiple sources of evidence, there ex- 
ist two models in DSmT [36,37]: 1) free combination model and 
2) hybrid combination model which is often used in real appli- 
cation because it takes into account some integrity constraints. 
In hybrid combination rule, it transfers partial conflicting beliefs 
to the corresponding intersected elements, but this increases the 
uncertainty of fusion results. The Proportional Conflict Redistribu- 
tion 1-6 (PCR1-6) rules [37,39,40] provide proper conflict redistri- 
bution ways, and they proportionally transfer conflicting masses to 
the involved elements. 

The difference of PCR1-6 rules mainly lies in the redistribution 
of conflicts, and PCR5 is considered as the most precise redistribu- 
tion way [37,39,40]. The combination of two independent sources 
of evidences by PCR5 rule is given as follows [37,39,40] 


m1@2(Xi) = »- 


Y,ZeG® and Y,Z40 
YNZ=X; 


m4(Y)-m(Z) (5) 


oy E& -mM2(Xj) 
xj¢6° and if my (Xj) +m2(Xj) 
XjNX j=0 (6) 
_m2(Xi)* m1 (Xj) ] ne 
T S| Xi eG and X; #9 
0 xX;j=0 


M1@2(Xi) + 


mMpcr5(Xj) = 





where G® can been seen as the power set 2°, the hyper-power 
set D® and the super-power set S®, if discernment of the fusion 
problem satisfies Shafer’s model, the hybrid DSm model, and the 
minimal refinement ©" of © respectively and where all denomi- 
nators are more than zero and the fraction is discarded when the 
denominator of it is zero [37,39,40]. 
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Nevertheless, PCR5 rule still has some disadvantages, such as, 
firstly, it is not associative in the fusion of multiple (more than 2) 
sources of evidences, so the combination order may have influence 
on the results, secondly, its computational complexity increases 
exponentially, when the focal elements’ number increases. Our re- 
search in this paper is mainly for reducing the complexity of PCR5 
within DSmT framework. 


2.3. The dissimilarity measure method of multi evidences 


The dissimilarity measure method of multi evidences and sev- 
eral Evidence Support Measure of Similarity (ESMS) functions have 
been given in [31,41]. The often used Euclidean ESMS function and 
Jousselme ESMS function are briefly recalled. 

1) Euclidean ESMS function Sime (m1, m2) 

Let © = {6;,62,---, On}, n > 1, |G®| be the cardinality of G®, 
m,(-) and m2(-) be two bba’s. The Euclidean ESMS function is de- 
fined by [31] 





|G?| 


Sime(m1,ma) =1— =, Ll —m2(Xi)]° (7) 


2) Jousselme ESMS function Sim ;(m;, mz) 
The Jousselme ESMS function [31] is defined based on the Jous- 
selme et al. measure [41] 





Sim j(m1,m2) = 1— + Vom — mz)" Dim; — m2) (8) 
J2 

where D = [Djj] is a |G°| x |G®| positively definite matrix, and 

Di =|Xi Xj|/[Xi UX; with X;, Xj ¢ G8. 

Some more ESMS functions can be seen in [31] for details. 
Simg(m1,mz) is considered with the fastest convergence speed 
[31], and it is adopted here as the dissimilarity measure for com- 
parison of the method proposed in this paper with the other meth- 
ods. 


3. An evidence clustering DSmT approximate reasoning method 


3.1. Mathematical analysis of PCR5 formula 


: . my (X;)?-m2(Xj;) 

As shown in Equation (6), D xj<c® ee ij Lima) + 
XjNXj=4 

mg (Xi)?-my (Xj) 

mig Xi) (X) | has symmetry. 


: sey M(Xi)?-ma(X;) 
Due to the symmetry, one item imi (Xi) m2 (Xj) 


Let m,(X;) =a and m2(X;) =x get 


is analyzed. 


m4(Xj) +m2(Xj)  a+x a+x 


Let X1,X2,°-°+,Xn € {m2(X ix j, Xj €G° and Xx; X; =}, then 


y E (Xj)? - md) 
m4 (Xj) +m2(Xj) 


XjeG® and i¢j 











XNY=6 
1 1 1 
=0)n o( + fovee op )} (10) 
a+x, ada+x2 a+Xn 
Let f(x) = a since f(x) is continuous function on (0, 1), it has 


a second order derivatives on (0,1), and f”(x) > 0 on (0,1), f(x) 
is a convex function. 

So 2(f (x1) + f (x2) +--+ fxn) = f RE), the equation 
holds iff xy = x2 =---=Xp. 


The approximate convex function formula is given by 


1 1 














a+x, a+ x2 a+ Xn 
n 
= +A, 
a+ (X1 +X2 +--+ +Xn)/n 
A>0, A=Oiffx; =x2 =---=Xp. (11) 








Let x1 < x2 <-+-<Xj <-+-<Xp, Carry out analysis of convex func- 
tion formula errors 








1 1 
-|a5 Steer | 





1 1 
+ |e 
A+X2 at (X1+X2 +--+ +Xn)/n 





| } : | (12) 
A+Xy, at (Xy+X2+-+++Xp)/n 


Analysis of the i item in Equation (12). 
Let (x1 +X2+--+-+Xn)/n=Xo, then 


1 1 1 1 
G+x; at (xy txot-:-+Xn)/n atx atx 
By Taylor expansion theorem 

1 1 
a+xji a+xo 

f’” (6) 

= f' (xo) (Xi (x) — Xo)” + 


5 € (Xj, Xo) OF (Xg, Xi), (14) 





(13) 








(xj — Xo)? +5", 





Xo) + J eo) 


then Equation (14) is transformed to 


1 1 1 n 














atx, a+x2 A+X, a+ (Xi +X2+---+Xy)/n 
= f'(xo)[(%1 — Xo) + (2 — X0) + +++ + Hn — X0)] 
J Ore, Xo)? + (x2 — Xo)? + +++ + (Xn — X0)?] 
+ Yoox — xo)’. (15) 
i=1 


Since f’(xo)[(x1 — Xo) + (x2 — X0) +--+ + (Xn — X0)] = 0, then 
1 1 1 n 




















atx, a+x2 A+Xn, a+ (Xi +X2+-+-+Xpy)/n 
f"(X0) = “ 
= Die — 20)? +) [0% — x0)”, (16) 
i=1 i=1 
where 
n 
>| 0(x; — x0)? 
i=1 
f"Xo) 
a: [ (x1 — x0)? + (X2 — Xo)? +--+ + (Xn — X0)?] 
ff" (81) f”’”" (82) 
+ a (x1 — xo)4 + 7 (x2 — xo) ++ 
f'"6 ) 
A! (Xn —Xo)* +++. 


Analysis of | f(x)|, m= 2, 3,---,00 


reo) =|( 1 y" =m 1 yr Hi) 
a+x a+x , 
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then 
(m—1) (m) 
we = a (x — X9)™—1) — een —) (x — Xo)™ 











(m—2) 
) |(x — x0)""1| 


1 ae ee zs 
aaulax) [= x0)" | 


(m—2) 
) |(x — x9)"~7| 


"gail 1 
~ (m—2)!\a+xo 








1 1 
~ (m—2)!\a+x 


















































1 1 
1 18 
«| aoa larm) sl (18) 
If x < Xo, |x — Xo] =X0 —X < Xo +4, then 
(1) (x (m) (x 
ERO) ig a YT ae ss fi @0) (, = xo)” i (19) 
(m— 1)! m! 
If m>2, x> xo, X<a+2xo, then 
: |(x—x0)| > 0 (20) 
. X — Xg)| > 0. 
m—1 a+ xo o 
So if x1 <a+ 2x9, i=1,2,---,n, 
(m1) (x 2 
Le —xo)™!] > |" (x —x0)™|,m>2, (21) 
(m—1)! m! 
namely, 
f" (Xo) - 0) 
5) (xj — X0)?| > (xj — X0)?| > 
(m) 
= 1), — Xo)" |. (22) 
m! 








Neglect the fourth order item errors and more order item errors. 
For the third order item is odd number item, for each x;, i= 1, 
2,°+:,n, Foo) (Xj — xo)? can be positive and negative. Then the 
sum of the third order items is much smaller than the sum of the 
second order items if xj <a+2x9,i=1,2,---,n. Neglect the third 
order item and more order item errors if xj <a+ 2x9, i=1,2, 























on. 
So, 
1 ; 1 1 n 
a+x, rte : "a+X_ A+ (xX +X2 +++: +Xn)/n 
2 
Xj —X 
oe xo)? 2 _ Die i= %0)" xj <a+2x9. (23) 
2(a — xXo0)3 
Then 
a “1 | x2 1 ...4 Xn 
atx, a+Xx~ " a+Xn 
| Xr bXa te + Xn 
a+ (X1 +X2 +--+ +Xn)/n 
” n n 2 
x «4 (Xj —X 
pal ( 0) > @i Xp)? = a2 eta o ; (24) 
ee 2(a+ Xo) 
From the above analysis, the errors are related to ee (Xj —X9)2 
az 
and 2(a+x0)3 


By the properties of if the mean point Xo increases, 


eae 
2(a+xo)>” 


2 . . 
Seca decreases quickly accordingly. When the cluster set {x;} 


is not particularly divergent, oy (xj — Xp)? is much smaller than 
divergent cluster. So get the conclusion that if the distribution of 
the cluster set {x;} is concentrated and the mean point Xp is large, 
the errors can be smaller. 

Based on the above error analysis, for reducing approximate 
error of the DSmT approximate convex function formula, a new 
Evidence Clustering method is proposed as follows: 


1) Force the mass assignments of focal elements in the evidence 
to two sets by the standard of = 


2) If x= 2, x; is forced to one set, denoted by eae and the sum 
of mass assignments for {x}} is denoted by S;, the number of 
points in {xi} is denoted by nz; otherwise, x; is forced to the 
other set, denoted by {x? }. 


3) If x € bea pick the focal element x; with the maximal value 


ae 1-S 
Ximax: If Ximax = a u 


4) Go on the step 3), until Xj max < 





, Xi is forced to one set ix } 


20— SL); ; the sum of mass as- 





signments for {x?} is denoted by Se. 


After the above cluster steps in the evidence, the mass assign- 
ments of focal elements are forced to 2 sets denoted by {xi} and 
{x?}. Compared to Pek the distribution of {x?} may be more con- 
centrated. Compared to {RP}; the main point of {xi} is large. So 
carrying out the above evidence clustering method in front of the 
approximate convex function formula can make errors of approxi- 
mate results much smaller. 

Let m1(X) =4, X1, X2,---,Xn € {M2(Xj)|Xj € Go,iFj and X;N 
Xj; =}, the approximate convex function formula of DSmT + 
PCR5 is given by 


3 eee 
Xj Xj 
XjeG® and i¢j ma) + maXyp) 
XNY=6 


2 a-n 

=a‘\n tA (25) 
a+ (X1 +X2+-+-+Xn)/n 

Finally, analysis of relationship between the approximate com- 


putation item and its errors item is as follows 
y [= (Xi)? « mid] 
m4 (Xj) +m2(X; 
Xj¢G® and i#j a nD) 
XNY=6 


tH gee oe nen) 
=a|n—a spared 
a+x, a+x2 a+ Xn 


Boral a-n 
a+ (X1 +X2+-+-+Xn)/n 
opie 11% — feo 
2[a + (x1 +X2 +--+ +Xn)/nP 
sal x1 +X. +++ +Xn 
a+ (X1 +X2 +--+ +Xn)/n 
ois [Xi — (X1 +2 +--+ +. Xn)/n? 
T 
2[a + (x1 +X2 +--+ +Xn)/nP 
After evidence clustering method, the influence of numerator 
X1 +X2+---+X, to the approximate computation item and the in- 
fluence of numerator a 4(X%i — (X1 + X2+- ++ +4+-Xn)/n)? to errors 
item is much smaller than the HUES of their denominators. So 
Nobo 
the approximate computation item a Tl is mainly 
2 fy Pi — (01 +24: -+Xn)/n]? 
Taree oe ee 
properties of convex function, all the errors items of focal ele- 
ments are negative. Based on the above analysis, the normalization 




















(26) 

















proportional to the error item a 
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method is applied as the final step of the proposed method for er- 
rors redistribution. 


3.2. The process of evidence clustering DSmT approximate reasoning 
method 


Based on the mathematical analysis of PCR5 in Section 3.1, 
the process of Evidence Clustering DSmT Approximate Reasoning 
Method is proposed as follows 


Definition 1. Assuming the existence of a cluster set {x}, the defi- 
nition of the total number of {x} is Num({x}) = number{x}, the sum 
of each point in {x} is Sum({x}) = }\{x} and the mean point of {x} 
is Mean({x}) = ee 

1) Carry out Clustering Evidence method proposed in Sec- 
tion 3.1. Force mass assignments of focal elements in each evi- 
dence to two cluster sets. Giving an example of two evidences, 
denoted by x = {xi}, y = {yi}, i= 1,2,---,n. Then force each 
evidence to two cluster sets, respectively denoted by x! = xt, 
C=(e} Us =x, y= (ih y= ih yUy ay. 

2) Un-normalized approximate fusion results are obtained by 
the approximate convex formula as follows 














aie x?Sum({yt}/ Yi) x?Sum({ys}) 
Pat" xit+Mean({y!}/yi) ° xi +Mean({y3}) 
yp Sum ({xt}/Xi) y?Sum({x$}) 





yi t+Mean({x!}/x;) 
yi € {yt}, x © {xh} 


yj+Mean({xs})’ 


x?Sum({y}}) 
x;-+Mean({y?}) 
__ yzSum({xt}/xi) 
" yi+Mean({x?}/xi) 
yi € {yt}, x1 © {xP} 
MCONVEXi = (27) 


x?Sum({y?}/Yi) 
x;-+Mean({y>}/yi) 
ypSum({x}}) 
yi-t+Mean({xs})’ 





Xi Yi + 








x.y, + esumnly) xpsum((y?}/Yi) 
ii xjt+Mean({y!}) * x;-+Mean({y>}/yi) 
y?Sum({x+}) y?Sum ({x? } /X;) 





" 'yi-+Mean((x!}) 
Vi € {yp}. xi € {XP} 
x?Sum({y!}/yi) x?Sum({y? }) 
xi+Mean({y;}/yi) * xi-+Mean({y}}) 
y2Sum({x*}) y2Sum({x> }/xi) 
yitMean({xt}) °  yj+Mean({x?}/xi)’ 
Vi € Ly? }, xi € {xt} 
3) Final approximate fusion results are obtained by the normal- 
ization step as follows 


meu = (1 - > 


X=XjN Yj 
XeG® 


yjt+Mean({x>}/x;) ’ 





Xi Vit 





MCONVEXi 


Xx) JS 28 
ms ) 4 Mconvexi oi 


3.3. Analysis of computation complexity 


If there are 2 evidence sources, assume that all possible unions 
of the elements in the discernment framework are empty for con- 
venient computation of complexity. All singleton focal elements 
and multiple intersection focal elements have mass assignments in 
hyper-power sets of 2 evidences, denoted by Ge? = {01, 62, +--+, On, 
Og M+ ++ Op, +++, OC NOGN-+-MOe}, {a, b,c, d,e} €[1,---,n], n denotes 
the number of singleton focal elements, c denotes the number of 
multiple focal elements. First, the computation complexity of 2 ev- 
idence sources fusion based on DSmT + PCRS is analyzed. Then, 
the computation complexity of the same 2 evidence sources fusion 


problem based on the method in this paper is also analyzed. Com- 
putation complexity comparison of two methods can be obtained 
from the analysis. 

Assume that the computation complexity of one time multi- 
plication is denoted by K, the computation complexity of one time 
addition is denoted by >, the computation complexity of one time 
division is denoted by w and the computation complexity of one 
time subtraction is denoted by B. The computation complexity of 
the processing procedure of 2 evidence sources fusion based on 
DSmT + PCRS, denoted by opsmr[n], is given as follows 


Opsmr+pcrs{n + c] =[K + (4K +2 + 4z)(n+c—1)]|(n+c) 
—x(2K+2y~+4z)+yu 
= (4n+ 4c — 3)(n+c)K 
+ (2n+2c—2)(n+c)w 
+ (4n+ 4c — 4)(n+c)X 
—x(2K+2y+4z)+ yu (29) 


where x denotes the number of multiple focal elements in the 
results, y denotes the number of the same multiple focal ele- 
ments generated in the procedure of mass assignments combina- 
tion product. 

The computation complexity of the same problem based on the 
method in this paper, denoted by ocy[n+ c], is given as follows 





oguin+c]=(n+c)K +2(n+0)[2GK+20+y)+ 2] 
+D4+ny+2(n+c)B—x(2K +2~+4zd)+ yz 
= 13(n+c)K + [4(n+c) + 1]= + (5n+ 4c) 
+2(n+c)B—x(2K +2W+4d)+ yd (30) 


where x denotes the number of multiple focal elements in the 
results, y denotes the number of the same multiple focal ele- 
ments generated in the procedure of mass assignments combina- 
tion product. 

Computation complexity comparison of two methods obtained 
from Equation (29) and Equation (30) shows that the computa- 
tion complexity of DSmT is almost proportion to (n +c)” and the 
computation complexity of the method in this paper is almost pro- 
portion to (n+). Analysis of computation complexity shows that 
the computational complexity of the proposed method increases 
linearly instead of exponentially when the elements number of the 
discernment framework increases linearly. 





4. Simulation experiments 
4.1. Simple cases of cluster sets in each evidence 


Example 1. If there are 2 evidence sources, assume that only 
singleton focal elements have mass assignments in hyper-power 
sets, denoted by Ge = {01,62,---,07}, k= or 2. The mass assign- 
ments in each evidence are x = {0.1, 0.1, 0.05, 0.3, 0.2, 0.2, 0.05}, 
y = {0.2, 0.05, 0.05, 0.2, 0.15, 0.3, 0.05}, the process of the method 
is given as follows: 

1) The mass assignments of focal elements in each evidence are 
clustered to two sets, denoted by x = {X1, X2, X3, X5, X6, X7} U {x4}, 
y ={¥1,¥2, 3,4, ¥5, 7} U {Yo}. 

2) Un-normalized approximate fusion results are calculated by 
Equation (27) 








MCONVEX1 = X1-° Y1 + y?  Surm( (x2, ¥3, 5» X6, *7)) Yi Xa 
yi + Mean({x2, x3, %5,X6,X7}) = V1 + X4 

x7 Sum({¥2, 93,4, V5, ¥7}) Xz Yo 

x1 + Mean({y2, ¥3,¥4,¥5,¥7}) x1+Ye 
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Table 1 
Euclidean similarities and average computing time comparison of the methods. 





Euclidean similarities with results of DSmT + PCR5 


Average computing 












































1 2 3 4 5 6 times (s) 
The method in this paper 0.9939 0.9911 0.9955 0.9938 0.9965 0.9960 0.0028 
The method in [32] 0.9583 0.9791 0.9588 0.9530 0.9443 0.9347 0.0105 
icine y3 - SuM({X1, X3, X5, X6, X7}) 5X4 In the same way, the Euclidean similarity between the results 
y3 + Mean({x1, X3, X5, X6, X7}) y3 +X4 MyXDL of the method in [32] and the results mpsmr+pcrs of DSmT + 
2 2 PCR5 is obtained by Equation (7) 
x5 -Sum({y1, ¥3, ¥4, ¥5, Y7}) x5 V6 
x5 + Mean({y1, ¥3. 4,5. ¥7}) %3+Ye Exp, = 0.9812. 
re - Sum({X1, X2, X5, X6, X7}) y3 -X4 From the above results of this example, the results obtained 
MCONVEX3 = 3° Y3 4 y2 + Mean({x;, X2,X5,X6.X7})  y2-+Xa by the method proposed in this paper have higher Euclidean Sim- 
3 fan 3 ilarity with DSmT + PCR5 than the existing DSmT approximate 
ne -Sum({¥1, ¥2, ¥4, Ys, Y7}) ne “V6 reasoning method in [32]. The Euclidean Similarity which remains 
2 + Mean({y1, y2, ya, Y5, Y7}) x2 + 6 over 99% shows that the method proposed in this paper has high 
5 accuracy and practical meaning. 
re a9) ged Yq: SuM({X}, X2, X3, X5, X6, X7}) 
eee ee y2 + Mean({Xx1, X2, X3,%5, X6, X7 Example 2. If there are the same 2 evidence sources in Exam- 
4 
2 2 ple 1, the hyper-power sets are denoted by G® = {6,,65,---, 6}, 
xq Sum(Y1, Y2, ¥3, 5, Y7}) X4° 6 k =1or2. The mass assignments in each evidence are x = 
x + Mean({y1, y2, ¥3, ¥5, Y7}) <1 + Ye {0.1, 0.1, 0.05, 0.3, 0.2,0.2,0.05}, y = {0.2,0.05, 0.05, 0.2, 0.15, 
y2 - Sum({x1, X2, X35 X6, X7}) y2-x4 0.3, 0.05}. The mass assignments in the second evidence source are 
MconvEXs = X55 + = See = unchanged, denoted by y = {0.2, 0.05, 0.05, 0.2, 0.15, 0.3, 0.05}, 
y5 + Mean({x1, X2,X3,%6,X7}) = Y5 +X and the mass assignments in the first evidence change and the 
x2 - Sum(Ly1, V2, ¥3, ¥4s Y7}) x %e sequence of the mass belief of each focal element moves one po- 
5 5 sition backward at one time to procedure 6 new evidences, such 
x5 + Mean({y1, Y2,¥3,¥4,¥7}) x5 + Ye as: 
ya - Sum({x1, X2, X3,X5, X7}) Ve Xa 
MCONVEX6 = X6- Y6 + —> a) x1 = {0.1, 0.1, 0.05, 0.3, 0.2, 0.2, 0.05}, 
ye t+ Mean({x1, X2,%3,%5,%7}) Ye +X 
2 X2 = {0.05, 0.1, 0.1, 0.05, 0.3, 0.2, 0.2}, 
xg Sum({y1, ¥2, ¥3, Ya, V5, Y7}) 
= {0.2, 0.05, 0.1, 0.1, 0.05, 0.3, 0.2 
x2 + Mean({y1, ¥2, ¥3. ¥4, Y5, Y7}) a a ae 
2 eae y2 - Sum({x1, X2, X3, Xs, Xe}) y2 as x4 = {0.2, 0.2, 0.05, 0.1, 0.1, 0.05, 0.3}, 
CONVEX7 = X7- Y7 
y3 + Mean({x1, x2, X3, Xs, X6}) y2 +%4 x5 = {0.3, 0.2, 0.2, 0.05, 0.1, 0.1, 0.05}, 
2 - Sum({y1, ¥2, V3, Y4 ¥5}) x V6 Xe = {0.05, 0.3, 0.2, 0.2, 0.05, 0.1, 0.1}. 





x5 + Mean({y1, y2. 3, ¥4,¥s}) %3+Ye 


Get Mconvex = {0.1588, 0.0558, 0.0273, 0.3108, 0.1926, 0.3108, 
0.0273}. 
3) The final approximate results are obtained by Equation (28) 


mcou = {0.1465, 0.0515, 0.0252, 0.2869, 0.1778, 0.2869, 
0.0252}. 


4) The results of DSmT + PCR5 and the method in [32] are also 
calculated for the comparison with the method in this paper. 
The results of DSmT + PCR5 


MpsmT+Pcrs = {0.1435, 0.0488, 0.0237, 0.2922, 0.1751, 0.2929, 
0.0237}. 
The results of the method in [32] 
Mxp_ = {0.1536, 0.0605, 0.0253, 0.2980, 0.1670, 0.2738, 
0.0217}. 


The Euclidean similarity between results mcy obtained by the 
method in this paper and the results mpsmr+pcrs Of DSmT + PCRS, 
is obtained by Equation (7) 


Ecu = 0.9932. 


Each new evidence and the evidence y are calculated to obtain 
the fusion results by DSmT + PCR5 and the approximate fusion 
results by the method in this paper and the method in [32]. Then 
Euclidean similarities of the approximate results of different meth- 
ods with the results of DSmT + PCR5 are obtained by Equation (7) 
and the average computing times of the methods are also taken 
record as Table 1. (In this paper, all the simulation experiments 
are implemented by Matlab simulation in the hardware condition 
of Pentimu(R) Dual-Core CPU E5300 2.6 GHz 2.59 GHz, memory 
1.99 GB.) 

As shown in Table 1, in the simple cases of cluster sets in each 
evidence, the accuracies of the method in this paper all remain 
over 99% and much higher than the method in [32]. Average com- 
puting time of the method in this paper is much lower than the 
method in [32]. At the same time, the accuracies of the method 
in this paper in different evidence cases change little, which prove 
that the method in this paper has higher performance stability. 


4.2. Complex cases of cluster sets in each evidence 


Example 3. If there are 2 evidence sources, assume that only sin- 
gleton focal elements have mass assignments in hyper-power sets, 
denoted by Ge = {01,62,---,012}, k = 1 or 2. The mass assign- 
ments in each evidence are x = {0.3, 0.35, 0.05, 0.05, 0.04, 0.06, 
0.02, 0.01, 0.02, 0.01, 0.04, 0.05}, y = {0.2,0.05, 0.04, 0.21, 0.15, 
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Table 2 
Euclidean similarities and average computing time comparison of the methods. 





Euclidean similarities with results of DSmT + PCR5 


Average computing 

















1 2 3 4 5 6 7 8 9 10 11 times (s) 

The method in this paper 0.9987 0.9983 0.9982 0.9979 0.9981 0.9985 0.9985 0.9983 0.9984 0.9983 0.9983 0.0038 

The method in [32] 0.8795 0.9330 09514 0.9484 08112 08636 0.8253 0.8342 0.8331 0.8189 0.8483 0.0186 
0.25, 0.05, 0.01, 0.01, 0.01, 0.01, 0.01}, the method process is given Table 3 
as follows: The mass assignments of highly conflict evidence sources. 

1) The mass assignments of focal elements in each evidence are Conflict evidence sources a b c d 

clustered to two sets, denoted by x = {X3, x4, X5, X6, X7, X8, X9, X10, Sy X—6 & 1—x-e € 
X11, X12} U {X%1, Xa}, VY = {¥2, 3, V7, V8, Yo. Yio, Y11, 12} U {¥1, 4, Sy € y-e & 1l-y-e 


y5, V6}. 
2) Un-normalized approximate fusion results are calculated by 
Equation (27) 


Mconvex = {0.3069, 0.2559, 0.0247, 0.1309, 0.0834, 0.1662, 
0.0177, 0.0019, 0.0041, 0.0019, 0.0110, 0.0153}. 
3) The final approximate results are obtained by Equation (28) 
mou = {0.3009, 0.2509, 0.0242, 0.1283, 0.0818, 0.1630, 0.0174, 
0.0019, 0.0041, 0.0019, 0.0108, 0.0150}. 


4) The results of DSmT + PCR5 and the method in [32] are also 
calculated for the comparison with the method in this paper. 
The results of DSmT + PCR5 


MpsmT+Pcrs = {0.3019, 0.2524, 0.0235, 0.1282, 0.0811, 0.1635, 
0.0169, 0.0018, 0.0039, 0.0018, 0.0104, 0.0146}. 
The results of the method in [32] 
mMxp. = {0.3710, 0.1834, 0.0276, 0.1269, 0.0828, 0.1651, 0.0003, 
0.0000, 0.0001, 0.0001, 0.0002, 0.0002}. 


The Euclidean similarity between the results mcoy and 
Mpsmr+pcrs is calculated, denoted by Ecy = 0.9984 and computing 
time is 0.0035 s. 

In the same way, the Euclidean similarity between the results 
Myxp_ and mMpsmt+pcrs is calculated, denoted by Exp, = 0.9287 and 
computing time is 0.0185 s. 

As shown in the above experiment results, the results obtained 
by the method proposed in this paper have higher Euclidean simi- 
larity and lower computation complexity than the existing method 
in [32]. The Euclidean similarity which remains over 99% shows 
that the method proposed in this paper has higher accuracy and 
practical meaning. 


Example 4. If there are the same 2 evidence sources in Example 3, 
the hyper-power sets are denoted by Ge = {01,02,---, 012}, k= 1 
or 2. The mass assignments in the second evidence source are 
unchanged, denoted by y = {0.2, 0.05, 0.04, 0.21, 0.15, 0.25, 0.05, 
0.01, 0.01, 0.01, 0.01,0.01}. The mass assignments in the first 
source change and the sequence of the mass beliefs of focal el- 
ements moves one position backward at one time to procedure 11 
new evidences. 

Each new evidence and the second evidence y are calculated 
to get the fusion results of DSmT + PCR5 and the approximate re- 
sults by the method in this paper and the method in [32]. Then 
Euclidean similarities of the approximate results of different meth- 
ods with the results of DSmT + PCR5 and the average computing 
times of the methods are taken record as Table 2. 

As shown in Table 2, under complex cases of cluster sets in each 
evidence, the accuracy of the method in this paper also remains 
over 99% and much higher than the method in [32]. Average com- 
puting time of the method in this paper is lower than the method 
in [32]. At the same time, the accuracies of the method in this pa- 





Euclidean similarity 








Fig. 1. Euclidean similarity of the method in [32] with DSmT + PCRS. 


per in different evidence cases change little, which prove that the 
method in this paper has higher performance stability. 


4.3. Cases of highly conflict evidence sources 


Example 5. In order to verify that information fusion of highly con- 
flict evidence sources can be effectively solved by the method in 
this paper, assume there are two highly conflict evidence sources 
with the hyper-power set denoted by D® = {a,b,c,d} and the 
mass assignments of two evidence sources are shown as Table 3. 

Let € =0.01, x, y € [0.02, 0.98]. The fusion results are obtained 
by different methods when x, y is increasing from 0.02 to 0.98 by 
0.01 step at the same time. Euclidean similarity of the method in 
[32] with DSmT + PCR5 is shown in Fig. 1. Euclidean similarity of 
the method in this paper with DSmT + PCRS5 is shown in Fig. 2. 

The average Euclidean similarity of the method in this paper 
is 0.9873 and the average Euclidean similarity of the method in 
[32] is 0.8513. It’s shown that information fusion problem of highly 
conflict evidence sources can be effectively solved by the method 
proposed in this paper. 


4.4. Convergence analysis 


Example 6. If there are 2 evidence sources, assume that only sin- 
gleton focal elements have mass assignments in hyper-power sets, 
denoted by Ce = {01,02,---,012}, k= 1 or 2. The mass assign- 
ments in each evidence are x = {0.1, 0.01, 0.02, 0.25, 0.15, 0.05, 
0.1, 0.1, 0.1, 0.05, 0.05, 0.02}, y = {0.5, 0.35, 0.02, 0.02, 0.01, 0.01, 
0.01, 0.01, 0.01, 0.01, 0.01, 0.04}. First, fusion results of two evi- 
dences x and y are obtained by different fusion methods. Then, 
fusion of the prior fusion results with y is carried out repeatedly. 
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Table 4 
Fusion results comparison of the methods for convergence analysis. 





Fusion results 





1 DSmT + PCR5 0.4287, 
The method in [32] 0.4067, 
The method in this paper 0.4285, 
2 DSmT + PCR5 0.5829, 
The method in [32] 0.5687, 
The method in this paper 0.5705, 


3 DSmT + PCR5 0.6439, 
The method in [32] 0.6171, 
The method in this paper 0.6118, 
4 DSmT + PCR5 0.6726, 
The method in [32] 0.6152, 
The method in this paper 0.6264, 


0.2600, 0.0058, 0.1163, 0.0539, 0.0109, 0.0297, 0.0297, 0.0297, 0.0109, 0.0109, 0.0136 
0.0994, 0.0212, 0.1888, 0.0798, 0.0385, 0.0027, 0.0027, 0.0027, 0.0010, 0.0010, 0.0019 
0.2610, 0.0060, 0.1151, 0.0535, 0.0110, 0.0296, 0.0296, 0.0296, 0.0110, 0.0110, 0.0141 


0.3324, 0.0031, 0.0388, 0.0120, 0.0017, 0.0053, 0.0053, 0.0053, 0.0017, 0.0017, 0.0099 
0.1686, 0.0036, 0.0642, 0.0202, 0.0127, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
0.3459, 0.0031, 0.0371, 0.0116, 0.0018, 0.0052, 0.0052, 0.0052, 0.0018, 0.0018, 0.0109 


0.3342, 0.0018, 0.0084, 0.0016, 0.0006, 0.0008, 0.0008, 0.0008, 0.0006, 0.0006, 0.0060 
0.1968, 0.0009, 0.0109, 0.0027, 0.0021, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
0.3655, 0.0021, 0.0078, 0.0015, 0.0006, 0.0008, 0.0008, 0.0008, 0.0006, 0.0006, 0.0069 


0.3181, 0.0012, 0.0017, 0.0004, 0.0004, 0.0004, 0.0004, 0.0004, 0.0004, 0.0004, 0.0040 
0.2047, 0.0008, 0.0018, 0.0004, 0.0003, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
0.3626, 0.0015, 0.0017, 0.0004, 0.0004, 0.0004, 0.0004, 0.0004, 0.0004, 0.0004, 0.0048 
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Fig. 2. Euclidean similarity of the method in this paper with DSmT + PCR5. 


Through the fusion results each time, the different methods’ con- 
vergences are analyzed in Table 4. 

As shown in Table 4, the convergences of three methods are 
similar. The results of the methods can converge to the main focal 
elements after 3 times of evidences fusion. However, the results 
of the method in this paper have higher Euclidean similarity with 
DSmT + PCRS, and lose less information than the method in [32] 
each time. 


4.5. Monte Carlo simulations in the case of non-empty multiple focal 
elements 


If there are 2 evidence sources, assume that singleton focal 
elements and multiple focal elements have mass assignments in 
hyper-power sets, denoted by P© = {6}, 02, +--+, 029,01 105 419M 
620}. Carry out 1000 times Monte Carlo simulation experiments. 
The simulation experiment process of each time consists of 3 steps. 
First, mass values of all focal elements of hyper-power sets in each 
evidence are assigned randomly. Then, the fusion results of 2 ev- 
idences are obtained by DSmT + PCR5 and the method in this 
paper separately. Thirdly, the Euclidean similarity of the method 
in this paper with DSmT + PCR5 and computing times of both 
methods are calculated and recorded. The Monte Carlo simulation 
results are shown in Fig. 3, Fig. 4 and Table 5. 

As shown in Fig. 3, Fig. 4 and Table 5, in the case of non- 
empty multiple focal elements, the average Euclidean similarity of 
the method in this paper can reach 98.49% and Euclidean similar- 
ity changes little with different evidences. Computing times of the 
method in this paper almost reduce halfly than DSmT + PCRS. 


-3 
s 
36 x 10 


—4— computing time of the method in this paper 
—8— computing time of DSmT+PCR5 
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Fig. 3. Computing time comparison of the method in this paper with DSmT + PCRS. 


4.005 —4:— Euclidean similarity of the method in this paper 
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Fig. 4. Euclidean similarity of the method in this paper with DSmT + PCR5. 


4.6. Monte Carlo simulations in the case of increasing focal elements 
number 


Example 7. If there are 2 evidence sources, assume that only sin- 
gleton focal elements have mass assignments in hyper-power sets, 
denoted by G® = {6),69,---, 19}, k= 1 or 2. The focal elements 
of the hyper-power sets are increased 10 focal elements each time 
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Table 5 
Fusion results comparison in the case of increasing focal elements number. 





Average Euclidean Max Euclidean 


Min Euclidean 


Average computing Max computing Min computing 




















similarity similarity similarity time (ms) time (ms) time (ms) 
DSmT + PCR5 1.9 31 1.9 
The method in this paper 0.9849 0.9956 0.9693 0.84911 1.4 0.83593 
4.005 —+.— Euclidean similarity of the method in this paper —4— computing time of DSmT+PCR5 
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Fig. 5. Euclidean similarity of the method in this paper with DSmT + PCR5. 


until 500 and the mass values of every focal elements of hyper- 
power sets in each evidence are assigned randomly each time. 
1000 times Monte Carlo simulation experiments in each hyper- 
power set are carried out and the average Euclidean similarities 
of the method in this paper with DSmT + PCRS are calculated and 
recorded as shown in Fig. 5. The average computing time com- 
parison of the method in this paper with DSmT + PCR5 in each 
hyper-power sets is shown in Fig. 6. 

The fusion results comparison in the case of increasing focal 
elements number is shown in Table 6. (As the increasing number 
of focal elements, the mass assignment of average cluster center is 
decreasing sharply. For reducing computation complexity, neglect 
the influence of the different classification of clusters and apply 
the standard of 2 as one step cluster method, n denotes the focal 
elements number of hyper-power sets.) 

As shown in Fig. 5, Fig. 6 and Table 6, 

1) In the case of increasing focal elements number, comput- 
ing time of the method in this paper decreases significantly, and 
the computation complexity of the method almost appears lin- 
ear growth instead of exponential growth, which proves that the 
method in this paper has a high application in the case of com- 
plex fusion problems with big number of focal elements. 

2) The accuracy of the method in this paper is increasing with 
the growth of focal elements number of hyper-power sets as the 


Table 6 
Fusion results comparison in the case of increasing focal elements number. 


the number of focal elements 


Fig. 6. Computing time comparison of the method in this paper with DSmT + PCRS. 


errors item becomes much smaller. The minimum average Eu- 
clidean Similarity is 0.9974 in the case of the minimum num- 
ber of hyper-power sets. When the number of hyper-power sets 
increases over 50, the average Euclidean similarity exceeds 99%, 
which proves that the method in this paper can effectively sup- 
port correct and quick decision in the case of large data. 


5. Conclusions 


A new Evidence Clustering DSmT Approximate Reasoning Me- 
thod is proposed based on convex function analysis in this paper. 
This new method consists of three steps. First, the belief masses 
of focal elements in each evidence are clustered by the Evidence 
Clustering method. Second, the un-normalized approximate fusion 
results are obtained using the DSmT approximate convex function 
formula, which is acquired based on the mathematical analysis 
of PCR5 rule in DSmT. Finally, the normalization step is applied. 
The method reduces computation complexity of PCR5 based in- 
formation fusion rule within DSmT framework significantly which 
constrains the wide application and development of DSmT and re- 
mains high accuracy. Simulation results in different cases show 
that the method in this paper can process evidences fusion prob- 
lems effectively and efficiently, especially, in the case of large data 





The focal elements number of hyper-power sets is increasing from 10 to 510 





Average computing times of 
DSmT + PCRS (s) 


0.0008 0.0032 0.0073 0.0132 0.0208 0.0297 0.0406 0.0530 0.0670 0.0829 0.1017 0.1219 0.1405 0.1636 0.1867 0.2133 0.2431 0.2685 
0.3009 0.3319 0.3698 0.4013 0.4427 0.4816 0.5212 0.5588 0.6011 0.6462 0.6959 0.7413 0.7909 0.8425 0.8962 0.9505 1.0082 1.0659 


1.1263 11882 1.2516 1.3166 1.3828 1.4509 1.5225 1.5926 1.6660 1.7409 1.8174 1.8946 1.9765 2.0574 2.1397 


Average computing times of the 
method in this paper (s) 


0.0003 0.0006 0.0010 0.0013 0.0017 0.0020 0.0024 0.0028 0.0032 0.0036 0.0040 0.0044 0.0048 0.0053 0.0057 0.0062 0.0067 0.0071 
0.0076 0.0081 0.0087 0.0092 0.0098 0.0103 0.0108 0.0114 0.0118 0.0124 0.0130 0.0136 0.0142 0.0148 0.0154 0.0160 0.0167 


0.0173 0.0180 0.0186 0.0193 0.0200 0.0207 0.0214 0.0222 0.0229 0.0236 0.0244 0.0251 0.0259 0.0265 0.0275 0.0282 


Average Euclidean similarities 
with results of DSmT + 


0.9774 0.9826 0.9871 0.9881 0.9895 0.9907 0.9917 0.9922 0.9922 0.9927 0.9931 0.9933 0.9939 0.9940 0.9938 0.9945 0.9945 
0.9946 0.9947 0.9947 0.9950 0.9953 0.9949 0.9954 0.9955 0.9955 0.9955 0.9956 0.9959 0.9958 0.9961 0.9959 0.9961 0.9960 


PCRS 0.9959 0.9961 0.9962 0.9964 0.9961 0.9964 0.9964 0.9964 0.9964 0.9965 0.9966 0.9966 0.9966 0.9968 0.9966 0.9968 0.9967 
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and complex fusion problems, the method can get highly accurate 
results and need low computation complexity. 

However, our proposed method is not associative and the com- 
bination order really has influence on the results, which is the 
same problem to PCR5. This problem is really needed to be re- 
solved especially for the efficient fusion of multiple (more than 2) 
sources of evidences and more research would be carried out to 
resolve this associative question of our method in the near future. 
Our method would also be applied to the realistic and applied data 
or database in future. 
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